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Abstract 

The collision operator for a relativistic plasma is reformulated in terms of an 
expansion in spherical harmonics. In this formulation the collision operator is ex¬ 
pressed in terms of five scalar potentials which are given by one-dimensional in¬ 
tegrals over the distribution function. This formulation is used to calculate the 
electrical conductivity of a uniform electron-ion plasma with infinitely massive 
ions. 
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I. Introduction 


Landau^ first obtained an accurate kinetic equation for a nonrelativistic plas¬ 
ma. The Landau collision operator was generalized to a relativistic plasma by 
Beliaev and Budker.^ The collision operator in both cases involves integrals of the 
distribution function of the background species. This makes the operators difficult 
to evaluate, numerically or analytically. 

In the nonrelativistic case this difficulty was removed by Rosenbluth, Mac¬ 
Donald, and Judd,^ and by Trubnikov.^ They recast the Landau operator into dif¬ 
ferential form by writing it in terms of derivatives of two scalar potentials. The 
potentials in turn satisfy a pair of elliptic partial differential equations. With the 
aid of this formulation, the numerical evaluation of the collision operator may be 
accomplished straightforwardly by solving the potential equations. An analytical 
solution of these equations in terms of spherical harmonics was given by Rosen¬ 
bluth et al.^ 

Recently,^ we formulated the relativistic collision operator of Beliaev and Bud- 
ker^ in terms of five scalar potentials, which again obey elliptic partial differential 
equations. In the present paper we extend this formulation by solving the potential 
equations in terms of an expansion in spherical harmonics. As an application, we 
calculate the electrical conductivity of a relativistic plasma with massive ions. 

In Sec. II we review the differential formulation of the collision operator. The 
spherical harmonic expansion is developed in Sec. III. It is shown in Sec. IV how 
the nonrelativistic results are recovered in the limit c ^ oo. We evaluate the 
zeroth-order and tirst-order spherical harmonic components of the collision term 
with respect to a Maxwellian background in Secs. V and VI. The calculation of the 
conductivity is given in Sec. VII. 

II. Potentials 

We begin by summarizing the differential formulation of the collision operator. 
This repeats our earlier exposition®; however, we now introduce a more compact 
notation and also write the results in such a way that the nonrelativistic limit is 
more easily obtained. 

The collision term for a plasma of species s colliding off species s' may be 
written in the Fokker-Planck form as^’^: 

fs') = ^ ^ j, (1) 
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where the diffusion and frietion eoeffieients are given by 


d*/"7u) = 


ps/s' 


F*/"'(u) = - 


2ns' 
pV^' rn 


\J{u,u')fs'{u')d^u, 

( d 


2ns' m,/ J \ du' 


■ U(u,u') ]fs'{u')(Pu\ 


( 2 a) 

(2b) 


ps/s' ^ Us'q^ql log 

dyrepm^ 

Here, u is the momentum per unit rest mass, and, in the relativistie ease, the kernel 
U is given by^ 


r2 

U(u, u') = —;—— uu — u'u' + r(uu' + u'u 


in whieh 7 = a/1 + u^/c^, 7 ' = a/1 + u'^/c^, and 

r = 77 — u • u /c , 


w 


= c\/r 2 - 1 . 


(3) 


(4a) 

(4b) 


The quantity r is the relativistie eorreetion faetor eorresponding to the relative ve- 
loeity of the two interaeting partieles. The relative speed of the interaeting partieles 
is given by m/r. In the nonrelativistie limit, r ^ 1 and w |u — u'|. 

In the previous paper^ we expressed the Eqs. (2) in differential form, making 
use of the potentials 


^s'[i]o(u) 

^s'[ 2 ] 02 (u) 

'^'s'[3]022(u) 

^s'[2]ll(u) 


1 

dvr j 

[w V.'K) , 


(5a) 

1 

Stt j 

[ y , 


(5b) 

1 

327r 

J sinh“^(m/c) - 

w c)fs'{u) , , 

) 7 ' 

(5c) 

1 

dvr j 

f -1 . / /A 

1 rw fs'{u) y , 


(5d) 

1 

Stt j 

1 csmh~^{w/c)fs'{u) 

So! 

i 

(5e) 
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These potentials satisfy the differential equations 


-^o'I^s'[i]o = fs', (6a) 

L2^s'[2]02 = 'l's'[l]0> (6b) 

-^2'I's'[ 3]022 = 'I's'[2]02) (6c) 

= fs'i (6d) 

-^i'I's'[ 2 ]ii = (6e) 


where 


La^ = I + 


uu 


3u 5^- 1-0^ 


5uclu ' c)u 


(V) 


In terms of these potentials the diffusion and friction coefficients are given by 


D*/"'(u) = - 


dyrT"/"' 


n,/ 


1 (, I uu\^ 

-(^'- + ^ + ^j'h.q2]02 

4 /, I uu\ 

-- L - ^ ^ U-.,, 


7c^ 


c2 “ I '*^s'[3]022 


F*/^'(u) = - 


TUs 1 
Us' mg' 7 


- ^K4'^q2]ii 


where 


, . f, uu\ 5^4' 

LT- u =1 + ^ 


(9uc)u 




T/ \ I, 

K>P(u)= 1 + ^ . —, 


(8a) 

(8b) 


Equations (6) and (8) constitute the differential formulation of the relativistic col¬ 
lision operator. Boundary conditions on the solutions of Eqs. (6) are obtained by 
expansion of Eqs. (5) for |u| ^ oo. 

The present notation differs from that in Ref. 5. The correspondence is 'hs'[i]o 
= ho, 4's'[2]02 = hi, 4's'[3]022 = h2, = go, 'I's'pjii = Pi- The present 

notation reflects more clearly the structure of the differential equations (6). By 
adjoining to Eqs. (5) the definition 


^s'[ 0 ](u) = /s'(u), 
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and by using in the context to stand for a string of A: > 0 indices, we 

can express the five potential equations (6) in the concise form 

La^s'[k+l]*a ( 9 ) 

The integral representations (5) will likewise condense into a single formula 

1 /* u.^ 

^s'[fc]*(u) = (10) 

for A: > 0. The kernel function yo[fc]* is defined in fhe nexf secfion. 

III. Solution to potential equations 

In fhis section, Eq. (9) will be solved by separafion of variables in spherical 
coordinafes. The choice of spherical coordinafes is a nafural one since frequenlly 
fhe disfribufion function is nearly spherically symmefric and so is well represenfed 
by only a few spherical harmonics. Since depends on fhe disfribufion of 

species s only, we will simplify fhe nofafion by dropping fhe species subscripfs. 


A. Spherical harmonic expansion 


In a spherical {u, 9, (j)) coordinafe sysfem, fhe operator La is 


= Y 


02 q, 

dv? 


+ 


3u\ d'if 1 / 

'2 ' (9ri u2 1 Q02 


+ cot 6 


u c 

1 


5T' 


+ 


92 ^ 1 - 
+- 


sm 


9 d(j)' 


Eel us expand fhe pofenfials in terms of Eegendre harmonics: 

00 I 

^[fc]*(u) = Xl Ym[k]*{u)Pr{cos9)e^p{im(l)). 

1=0 m=—l 


( 11 ) 


( 12 ) 


The coefficients are given by 

(1 — 777 V /* d^\ 

Ym[k]*{Y) = (21 + 1 ) ^^ J —^[k]*{Y)PY{cos9)ex.p{-im<j)), 

where dQ is an element of solid angle 
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Equation (9) becomes 
where 


^l,a'4^lm[k+l]*a 


Ll,aX{u) 


1 + 


u 


d'^X 

du^ 


+ - + 
\ u 


3u\ dx 


l{l + 1 ) a? — 1 


X, 


(13) 


(14) 


and where '4^ira[o] (u) = fim{u) is the {I, m) coefficient in the expansion of /(u) in 
Legendre harmonics. 

As before, Eq. (13) stands for a set of five differential equations for fhe po¬ 
tentials V'zrnipo, V’«m[ 2 ] 02 > V’«m[ 3 ] 022 , V'Zmli]!, and We will, howcvcr, find if 

convenienf fo solve for fhe pofenfials wifh arbifrary indices, i.e., fo solve fhe sysfem 
of equations 


d-‘l,a'4^lm[l\a 

V’im[2]aa' 

Ll,a"'4^lm[3]aa'a" 


firm 

(15a) 

4lm[l]a7 

(15b) 

4lm[2]aa' > 

(15c) 


for arbifrary a, a', and a". 

Eor lafer reference we lisf fhe componenfs of fhe operafors L and K fhaf occur 
in Eqs. (8): 


L 'h 

^UU ^ 








^e4>^ 


4 0 ^ 4 ' 7 ^u O'f 

du^ du ’’ 

1 0^4' 7 ^ 04' 

V? 86“^ u On ’ 


1 


024- 72 1 

H--1-^ cot I 


V? sin2 9 d(jr u du u 
72 I Q‘2\^ I 

u 89 






= 


^8u89 
/ 0^4' 


7' 


04- 


1 04'\ 
u sin 9 8u84> u 8<p J’ 

1 / 0^4^ 04^^ 

n2 sin 9 y 8984) 94’) ’ 


(16a) 

(16b) 

(16c) 

(16d) 

(16e) 

(16f) 
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= 72 

Ke^' = - 

= 


9^- 

du ’ 

1 

u do ’ 

1 5^- 


iisin^ dcf) ’ 



(16g) 

(16h) 

(16i) 

(16j) 


B. Homogeneous solutions 

In order to solve the inhomogeneous equations (15), it is required first to deter¬ 
mine the solutions to the homogeneous equations 



(17a) 


(17b) 

^l,a"'^lm[S\aa'a" '^lm[2]aa' ‘ 

(17c) 


It is shown in the appendix that two independent solutions to Eq. (17a) are 

V’L®[i]a(«) = 7qi]a(«/c) and V’L®[i]a(^) = yma{u/c) where 


ma{u/c) = y (18a) 

yi[l]a{u/c) = (-l)-^-^y^P^+J/2(7)> (18b) 

and Pjf is the assoeiated Legendre funetion of the first kind. 

When the full system of equations (17) is eonsidered, we need to introduee the 
additional funetions, defined reeursively by Eq. (A 6 ), 


jl{k+2\*aa' 


j)[A:+l]>i!a 7)[A:+l]>i!a' 

a? — a'2 
dji[k+i] 
d{a?) ’ 


for a / a', 
for a = a'. 


(19) 
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with yi\k]* defined in a similar fashion. The general solution to Eqs. (17) is given 


1 

si 

1 _ 


1 

3" 


1- 

^lm[2]aa' 

II 

^ jl[2]aa' 


^ yi[2]aa' 

,/,HS 

l^lm[3]aa'a" i 


4 • 

^ Jl[3]aa'a'' 

_C^yi[3]aa'a''_ 



1 

O 

1_ 


1- 

o 

1_ 

+ Qm[2]aa' 

jl[l]a' 

^lm[2]aa' 

yi[l]a' 


_C^jl[2]a'a" _ 


_C^yi[2]a'a'' _ 


+ Ci 


lm[S\aa'a" 


-1 

O 

_1 


1- 

o 

1_ 

0 

^lm[3]aa'a" 

0 





( 20 ) 


where and are arbitrary eonstants (independent of u) and the argu¬ 

ment of the funetions jiyk]* and is u/c. 

The funetions and are invariant under permutation of the indiees in 
* and invariant under ehange of sign of any index in Also the funetions satisfy 

ji[k]*{-z) = {-lyji\k\*{z) and yi\k]*{-z) = {-Vf yi\k\*{z)\ and and 

are related by 


It follows from Eqs. (All) that c^^^^~‘^ji[k]*{u/c) and reduee 

to finite and nonzero expressions in the nonrelativistie limit: 

,,/+ 2 fe -2 

lim = 


u 


lim c 

C^OO 


(2A:-2)!!(2/ + 2A:-l)!!^ 
.2k-i-3 / / ^ ^ (-l)^(2t - 2fc + 1)!! 

yi[k]*[u/c) (2A:-2)!!u'-2*:+3 ' 


(21a) 

(21b) 


Eurther properties of these funetions are given in the appendix. 

Eor our problem, I and the indiees * are always integers, in whieh ease 
and yi[k]* may be expressed in terms of elementary funetions. The funetions that 
we need are given explieitly in Eqs. (A27) and (A28). It is seen there that the 
kernels appearing in Eqs. (5) are preeisely the funetions yo[k]*', this justifies the 
general definition for the potentials given in Eq. (10). 


C. Green’s function 

With the solutions of the homogeneous equations in hand, it is a straightforward 
matter to eonstruet a Green’s funetion and so to write down the general solution of 



the inhomogeneous problem. The final task will be to apply appropriate boundary 
eonditions from Eq. (10). 

We begin by defining the funetions 


) = 0, (22a) 

Nl[l]a{u,u') = c~^yi[l]a{u/c)jili]a{u'/c), (22b) 

Nl[ 2 ]aa’{u,u) = c(yi[l]a{u/c)ji[ 2 ]aa’{u / c) + yi[ 2 ]aa'{u/c)ji[i]a'{u / 

Nl[ 3 ]aa'a"(.U,u') = (y^lja (^^/c)jz[3]aa'a" + yq2]aa'(^^/c)j)[2]a'a" («7c) 

+ yi[3]aa'a"{u/C )ji[i]a"iu'/c)). (22d) 

From Eqs. (21), we see that the funetions Ni^k]* reduee to finite nonzero expres¬ 
sions in the nonrelativistie limit. These funetions satisfy the differential equations 

L'l,a^l[k+l]*a(y^ )'^) ^l[k]*i'^ 

Next we eombine u') and u) to obtain the Green’s funetions 


where n> 


Ki[k]*{u,u') = 
max(n, u') and tt< 


^6(u — u'), for A: = 0 

(23) 

Ni[k]*iuy,u<), fork>0, 

min(n, u'). These Green’s funetions satisfy 


To establish this relation, we use 


A 

du 


{Ni^k]*{u,u') - Niyk]*{u\u)) 

A 


II 


1 

'yu 

0 , 


for k = 1, 
otherwise. 


where the result for A = 1 follows from the expression for the Wronskian, Eq. 
(A5). Furthermore, the Green’s funetions satisfy equations analogous to Eqs. (19), 
namely. 


K, 


l[k+2]^aa' 


a? — a'2 
d{a?) ’ 


for a / a', 
for a = a'. 


(24) 
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The Green’s funetions therefore symmetrie under interehange of the in- 

diees in 

A partieular integral for the differential equations (13) is now given by 



We shall show next that this partieular integral is preeisely the solution defined 
by the integral form of the potential Eq. (10). This will be done by matehing the 
behavior of Eq. (25) near n = 0 to the behavior of the solutions of Eq. (10). Eor 
tt ^ 0 we have tt') ^ u). By using Eqs. (22) and (Alla) we 

obtain in this limit 

l 2k—1—3 pco 

^ ^21 + 1 )!! Jo ( 26 ) 


The general solution is obtained by adding the homogeneous solution, Eq. (20), 
to the partieular integral: V'/rnffc]* = It remains to determine the 

eoeffieients and appearing in Eq. (20). The primed eoeffieients 

C'irn[k]* found by eonsidering just the order of growth of 'ipim[k]*{u) at u = 0. 
Erom Eq. (12) it follows that = 0{v}), whereas from Eqs. (20) and 

(All), we have The eoeffieients must vanish 

in order to suppress this divergenee at u = 0. To establish that the unprimed 
eoeffieients also vanish, we must eonsider more earefully the behavior of 

i’im[k]*{u) at the origin. 

Erom Eqs. (20), (26), and (Alla), we find that the leading order behavior of 

'4^lm[k]* 


'^lm[k]* 


Clm[k]* T 0 


2k-3 



{u/cY 
(2( + l)!!' 


(27) 


In order to determine we expand the integral representation of riear 

u = 0. Substituting the spherieal harmonie expansion, Eq. (12), into Eq. (10), we 
obtain 




= 




flmju') Yim{6'A') 
l' Yim{0,4>) ’ 


(28) 


where we have introdueed the spherieal harmonies® 


Ylm{e,(t2) 


21 + I {I - m)\ 
dvr {I + m)\ ^ 


(cos 6) exp{im(p). 
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Because we want to expand Eq. (28) for small u, we write w from Eq. (4) as 


w = 



2eu'y, 


where 


e = rtcos a — (7 — 1 )^, 

v! 


and a = cos“^(u • u'/rttt') is the angle between u and u'. Because the u de¬ 
pendence only enters through e and because e = 0{u), we proceed by expanding 
Eq. (28) for small e. We use Eq. (A24) to expand yo[fc]*(^/c) as 


00 i' 

-i' ^ ^ e 


yo[k]*{w/c) = yi'[k]*{u'/c)jjj 

l'=0 ^ ■ 


(29) 


When we substitute for e and perform the angle integrations, we encounter integrals 
of the form 

Ini = COS^a 

Using formula (7.126.1) of Gradshteyn and Ryzhik,^ we can expand cos"a as 


cos’^a = 


n! 2n — 4A: -I- 1 
2>^kl + 


Eikewise, we can use Eq. (3.62) of Jackson® to expand P„(cos a) as 


Pnicos a) 


47r 

^ Y:,{9'A')Ynk{e,<p). 


Substituting these series into I^i and using the orthogonality condition for the 
spherical harmonics 


ldnYir^,{9,(j))Yimi9A)=Sin6„ 


we obtain 


^nl 


n! 

2 (n- 0 / 2 ((n-))/ 2 )!(n + / + l)!!’ 

0 , 


forn > / and n 
otherwise. 


I even, 
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Therefore, the first term in the sum in Eq. (29) which contributes to the integral in 
Eq. (28) is the term I' = 1. This results in a term of order uK The remaining terms 
in the sum contribute terms of higher order in u. Thus, we have 

yl ^2k—l—3 POO yj2 

V’«m[fc]*(«) ^Vo yi[k]*{u'lc)—fim{u')du'. (30) 


Comparing Eqs. (27) and (30), we find fhaf Qmffc]* = 0, and fherefore fhaf 
~ desired solution is given jusf by fhe particular integral, Eq. (25), 

which for /c > 0 we may wrife as 



wifh given by Eqs. (22). This completes fhe solufion for fhe pofenfials 

In obfaining fhis resulf, we have, in effecl, found a spherical harmonic decom- 
posifion of fhe kernel in Eq. (10) 


OO I 

;=0 m=—l 

OO 

= ^(2f + l)fV^fc]*(w>,u<)P/(cosa). 
1=0 


Eor k = 1, fhis yields an addifion formula for a class of associated Eegendre 
functions: 


(^2 _ 1 ) 1/4 


— ^(“ 1 )^( 2 ^ + 1)\/17 


1=0 


K-i/Y) 


K-uYy') 


2 (72 - 1)1/4 (y2 _ 1)1/4 


Pi (cos a ), 


where r = 77 ' — 1/72 — 11 / 7^2 — 1 cos a and 1 < 7 ^ < 7 . This identify is nof 
found in fhe usual handbooks, allhough if did furn oul lo be known.®’® 


IV. Nonrelativistic limit 

Taking fhe limil c ^ 00 in fhe preceding equations, we recover fhe well-known 
nonrelafivisfic collision operator. We cafalog here fhe imporfanf resulfs. 
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The kernel U reduees to the one given by Landau^ 

U(u, u') = 1 - — u'p I — (u — u')(u — u')^. 

The quantities r and w reduee to 1 and |u — u'|, respeetively. The potentials 'I'sqfc]* 
only depend on the number of indices, k, and not on their values. We will, there¬ 
fore, drop the indices and write instead of When Eq. (21b) is 

substituted into Eq. (10), the potentials become 




fs'{u), fork = 0, 

1 /* 111 _ 11 ^ 12^—3 

-Jil {2k-2)< /-'("W. fort>0. 


The potentials satisfy where L is now the velocity-space Eapla- 

cian, 

^ ^ d 

■ -7^- 

du du 


In particular, we have 


= fs', 

L'^s'[2] = 


The diffusion and friction coefficients are given by 

^ dyrT"/"' 92 

D*/" u =- 

Us' ouou ^ ' 

p,/.',,,. 4irr”-'*' 7n, d 


The homogeneous solutions to the separated radial components of the nonrela- 
tivistic potential equations are given by Eqs. (21). Substituting these into Eqs. (22) 
and then substituting the result into Eq. (31), we obtain the Eegendre harmonic 
expansion for the potentials 


^ ^ ^ 12 . f t l\ ! 


1 


u 


21 + 1 


du', 
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u'‘^flm{u') du' 


'>Plm[2] (u) = 




1 

1 

1 

u'^ 

2(2/ + 1) 

2/ + 3 

2/ - 1 u'-i 

1 

1 

1 

u^ 

2(2/ + 1) 

2/ + 3 

2/ - 1 u''- 


u'^ flra{u') du'. 


These expressions coincide with the results of Rosenbluth, MacDonald, and Judd.^ 


V. Isotropic background 

In this section and in Sec. VI, we consider the cases where the background is 
described by the I = 0, m = 0 and / = 1, m = 0 components in Eq. (12). First let 
us consider the azimuthally symmetric case, m = 0 and I arbitrary. For simplicity 
we will drop the m subscript and thus write 

OO 

^[fc]*(u) = J2^m*iu)Pi{cose). 

1=0 

When we substitute this representation into Eqs. (16), the second derivatives in the 
component may be eliminated by using the differential equation (13); this 

gives 


i(V’qfc+i]*a(«)^z(cos 6 »)) = 

'/(/ + !) a^-1 


2 , 27 ^ di^iik+i], 

7 i^i[k]* - 

^ ' u ou 


+ 7 




i’llk+i]* 


Piicose). (32) 


Fet us now specialize to an isotropic background I = 0, Po(cos 9) = 1. If we 
substitute Eqs. (16) and (32) into Eqs. (8), we obtain 


s/s’ ^ 47rr^/^' 7 

““’0 n,/ u 


2 dll^s'0[2]02 , 

27 - ^ -^^V’s'0[l]0 


87 ^ d'lps'0[3]022 8u 

du + ^^*'0[3]022 , 


(33a) 
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s/s' _ 47rr^/^' 1 

Us' -JU 


-7 


2 d'lps'0[2]02 U 


-2'*^^'0[2l02 


47^ d'lps'0[3]022 4tt 
+ ^^s'0[3]022 


_ 47rr®/^' rus 

^u,0 — 


n,/ m,/ 


■7 


'^V's'0[l]l _|_ 2 dV's'Op]!! 


dtt 




(33b) 

(33c) 


s/s' sis' 

The other components of Dg' and Fg' vanish. Lastly we substitute for the po¬ 
tentials using Eq. (31). Some massaging of the result leads to 


s/s' ^ 47rr*/"' [ rro.,2„2./ o„2./ 1 7 y: 

^uu,0 


n,' 


/2 


27 C Jg[ 2]02 8 c jo[ 3]022 3 / /s'o(^* ) 

/O'- u j 


+ 


27 '^c^J 0 [ 2]02 - 8 c 2 io[ 3]022 du \, (34a) 


-^06,0 ~ 


ns' Jo 


1 


1 


>-^ 0 [l ]2 „2 + ^2 P 0 [ 2]02 


7^ 


4 

~l” o o J 


^2 ^2-^0[3]022 


,'2 


-^fs'o{u')du' 
u 7 


+ 


1 


1 7 ^ 

2:^^ow 2-^1 ^ + :^)-?o[2]o2 


4 . 

+ ::^:j^-^o[3]o22 


u 


l—Js'o{u')du' 

7 


' 47rr^'^^' fTbs \ /■“ r 2 ■/ 




1 u‘ 


12 


(34b) 


Us! m-s' \ Jo L 

roo 


7 Jo[i]i 2 jg[ 2 ]iiJ ^2 y /s'o(n ) du 


+ / 4—io[2]02/s'o(nO <^n' 

./li W 


(34c) 


where jqfc]., = jim^u/c) and = ji[k]*{u'/c). 

Of particular interest is the case of a Maxwellian background, i.e., 


/s'o(n) — fs'm — 


Usinis' 


dvrcTs' K2 (m-s' IT s') 


exp 


ms'C ^ 
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where Kn is the nth order Bessel funetion of the seeond kind. First of all we ean 
verify that 

jps/s' _ TUsV j-^s/s' 

^ U.Q rp 

J^s' 

where v = n/ 7 . This is aeeomplished by substituting 

fs'm{u') = - ^ , fs'm{u') 

nis'v' aw 

sis' s/s' 

into the expression for q and integrating by parts. This relation between FJq 

s/s' 

and q implies that the collisions will cause fs to relax to a Maxwellian with 
temperature T^/. 

In the high-energy limit ms'C^{'y — 1) 3 > T^/, the indefinite limits in the inte¬ 
grals in Eqs. (34) can be replaced by 00 . We can perform the resulting integrals 
using formula (7.141.5) of Gradshteyn and Ryzhik^ which gives, after a change of 
integration variable, 

For I = 0 and a = 2 this gives the normalization condition for the Maxwellian 


roo 

/ A'KU^fs'm{u)du = ns'. 
7o 


On carrying out the integrations in Eqs. (34), we obtain 

j^s/s' _ ps/s' -^1 (I _ -^0 \ f 35 al 

^3X2!, iFl 72 c 2 j’ 

\ , KquI, 

^66.0 ~ ^ n ^ 2 ' 7 2 ]' 2 22 ’ 

’ 2 v K2 \ J K2 7 

= (350 

’ nis'K2\ ivi7^c^y 

where wfg = Tg/mg and the argument for the Bessel functions is rngW? jTgi. In the 
limit rrig' oo, we recover the Eorentz collision operator 

Duu,o = Fufl = 0) (36a) 

Dee,o = (36b) 

2 u 
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VI. First harmonic 


In addition to the isotropic components of the collision operator calculated in 
the previous section, the first harmonic of the Legendre expansion is also required 
in the calculation of the electrical conductivity. Specifically, we need fo compufe 
fhe ferm {fsm, fs'i cos 9 ). We can express fhis in terms of the potentials and 
their derivatives using Eqs. (1), (8), and (16) to give 


C^/^'{fsmiu),fs'liu)cOs9) 

fsm{u) COS 9 


4vrr"/"' 


Us' 


ms 

ms' 



u #5'1[1]1 
Uts du 



2 u d'lps'i[2]ii 
c^ufs du 


U #s'l[l]0 f u'^ 1 V 

Uts du ul) 


+ 



2 u \ dlps>l[2]02 
c^ul) du 



c^ul') 


'0s'l[2]O2 


8-fU dlps'l[3]022 
c'^ufs du 


+ 


87 

;CT^^'1[3]022 


(37) 


Finally, we substitute for the potentials and simplify to obtain 

C‘'/^'{fsm{u),fs’l{u) COS 9) 
fsmiu) COS 9 

dvrr^/^' f m, 1 . , 

= - < - fs'l{u) 

Us' ms' 7 


+ 


/o 


]_ f 2 ^ 7i[ 1]2 _ ^Q 7q2]02 


ms' 


u 


'ts 


U- 


ts 


I 1 ( ^ rns 7i[i]i ^ ^ ms 7i[2]ii 




m,/ W- 


ts 


ms' w 


ts 


+ 6 


*^^71 [2] 02 
Uts 


- 24 


'^^7i[3]022 
Uts . 
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where, as before, = jiyk]*{u/c) and = jiyk]*{u'/c). 

Both in Eqs. (34) and in Eq. (38), one ean substitute for jiyk]* from Eqs. (A27) 
and (A28) and thereby express the eollision operator entirely in terms of elementary 
funetions. The resulting expressions would be very badly behaved numerieally 
near u = 0 beeause of large eaneellations. It is preferable, therefore, to evaluate 
jiyk]* direetly by the method outlined in the appendix. 

The eollision operator obeys the eonservation law 

J [hsC^/^'ifs, fs') + hs'C^'/^fs'Js)] = 0, 

where hs = ao + sli ■ m^u + a 2 msc‘^'y, and oq, ai, and 02 are arbitrary eonstants. 
The eollision operator is also self-adjoint: 

j (xfsm, fs’ra) = J XC^ ^ {i’fsm, fs’ra) d^U, 

and satisfies the symmetry 

J 'll>C^^^'ifsm,Xfs'm)d^U = J xC'‘'^%fs'm,'ll’fsm)d^U, 

where V’ and x are arbitrary funetions of u, and = T^/. Combining these two 
properties gives 

C^^^'ihsfsm, fs'm) + hs^fs'm) = 0. 


This provides a useful eheek on the implementations of Eqs. (34) and (38) 


18 



VII. Calculation of the conductivity 


At this point the calculation of the electrical conductivity is straightforward. 
We consider an electron-ion plasma with infinitely massive stationary ions; m* ^ 
oo and fim ni5{u). In the presence of a weak electric field Ez, the electron 
distribution is given to first order by fem{^ + f) cos 9). The linearized Boltz¬ 
mann equation may be written in the form 


dxi 

dt 


QeEv 

Te 


+ Ce(Xl)) 


(39) 


where Ce{xi) is the linearized electron collision term 
1 


CeiXl) = 


fem COS 9 L 


C^^^UemXl COS 0 , fern) + femXl COS 9) 

+ C^I\femX\ COS 0, fim)] ■ 


The first term here is given by 


C ^ (/emXlCOS0 ) fern ) 1 ^ 2r.e/e , TT^eledXl 2 ^e/e 

fem COS 9 vfdu'' ^ du Vf 

with D and F given by Eqs. (34). The second is given directly by Eq. (38) with 
s = s' = e and fg'i = fs'mXi- The last term is given by the Eorentz limit, 
Eqs. (36), 

C^^KfemXlCOs9,fim) _ 

fem COS 9 

The conductivity is defined by 


a = 


47rge 

3F 


fem{u)Xliu,t 00)VU^ du. 


The time asympfofic solution f ^ oo) is defermined by solving Eq. (39) as 
an initial value problem. We lake Xi(^)^ = 0) = 0 (for example). The differ¬ 
ential terms ifemXi cos 9, fem) and {femXi cos 9, fim) are bofh freafed 
fully implicilly, while fhe integral term C^^^{fem, femXi cos9) is freafed explic- 
illy. This permils large time sfeps fo be faken and leads fo a rapid convergence fo 
a steady slate. A check on xi is oblained by evalualing fhe firsl momenl of fhe 
linearized Bollzmann equalion. The eleclron-eleclron collision terms drop out by 
conservation of momentum, to give 


^ [ femXll du 
o Jo 


rieqeE 
rUeZV^I^ ’ 


(40) 
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where 


Z = 


pe/i 

pe/e 


Qi 


log A®/* 


log A®/® 


<le 


is the effeetive ion eharge state, and where we have assumed qeUe + qiUi = 0. 
It is eonvenient to write a as 


m, 




A-Kel 


-. 3/2 


my^g^logA®/® Z 


-d(0,Z), 


where 

T T 

0 = —^ ^— 

nieC^ 511 keV’ 

The normalized eonduetivity d is a dimensionless funetion of two dimensionless 
arguments. In the limit 0 ^ 0, d is bounded and nonzero, and we reeover the 
nonrelativistie sealing a (x Te . Values of a for various 0 and Z are tabulated 
in Table 1 and plotted in Fig. 1. The nonrelativistie eonduetivity was first ealeu- 
lated by Spitzer and Harm,^'^ who quote values of d(0, Z)/a{0, oo). Their results 
eoineide with ours in the limit 0^0. 

In the limit Z ^ oo, eleetron-eleetron eollisions ean be ignored, and the rele¬ 
vant eollision term is the Lorentz eleetron-ion eollision term, Eqs. (36). This ease 
is eonsidered by Lifshitz and Pitaevskii.^^ We ean write 



At 


QeE 2 2 
^p^pe/e^ ^ 


The resulting eonduetivity is 


Att 1 
3 neufg 



U^V^fem du 


1 

307/2iT2(0-i) 



(f - 1 )^ 


exp 



dj. 


Evaluating the integral, we obtain 
1 

307/2^2(0-^) V 0 


- (1 - 0 + 202 


60^ - 240^ 


240^)exp(0 


-E 


(41) 
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where is the exponential integral. In the limit 0^0, this reduees to d = 

161 / 27 ^. For 0 ^ 00, we obtain a = 4/\/0, whieh agrees with the result of 
Lifshitz and Pitaevskii. 

Another traetable, albeit less interesting, limit is Z ^ 0. In this ease, this elee- 
trons equilibrate with themselves so that their distribution is a Maxwellian drifting 
at Vd and xi = Vdu/ute- The drift speed Vd is found by applying Eq. (40) to give 

^ _ 3exp(0~^)Ar2(0~^) qeEuj^ 

~ \/0(l + 20 + 202) meZr^/e ■ 

The resulting eonduetivity is 

^ 3exp(0-i)A:2(0-i) 

(J — —^-. (4^) 

\/0(l + 20 + 202) 

In this expression we reeognize the result obtained by van Erkelens and van Eeeu- 
wen^2 on the basis of a lowest-order variational treatment of the relativistie Boltz¬ 
mann equation. Their result for the eonduetivity of a relativistie plasma, therefore, 
eorresponds to the limit Z ^ 0. Eor this ease the limit 0^0 gives a = 
and the limit 0 ^ 00 gives d = 

VIII. Conclusions 

In our earlier work,^ we gave a differential formulation for the eollision op¬ 
erator for a relativistie plasma. This formulation is summarized by Eqs. (6) and 
(8). A major objeetive of the present work is to solve the differential equations (6) 
and, henee, to express the potentials in terms of quadrature. This was aehieved by 
using an expansion in Eegendre harmonies, Eq. (12). The radial eomponents of 
the potentials are then given by Eq. (31) where the kernels u') are given 

by Eqs. (22); these in turn involve the speeial funetions and whose 
properties are given in the appendix. The entire formulation is well-behaved in the 
nonrelativistie limit; indeed, in this limit, the potentials and their solutions agree 
with the earlier nonrelativistie treatment of Rosenbluth et al.^ and Trubnikov.^ 

Several eomputer eodes exist whieh solve the nonlinear Eokker-Planek equa¬ 
tion in the nonrelativistie limit. In many of these eodes the eollision operator is 
evaluated in terms of a Eegendre harmonie expansion of the potentials. Our results 
are easily ineorporated into sueh eodes, allowing them to treat relativistic colli¬ 
sions. The evaluation of the collision operator will be a few times more costly than 
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in the nonrelativistic case, firstly because five pofenfials need fo be compufed in- 
sfead of fwo, and secondly because fhe kernels involve fhe special funcfions 
and Hiyk]* insfead of simple powers of u. 

As an application of Ibis formulalion, we give in Eqs. (34) explicif forms for 
fhe diffusion and friction coefficienfs for an isofropic background. Finally, we 
calculate fhe elecfrical conducfivily of an elecfron-ion plasma wifh massive ions. 
Our resulfs agree wifh fhose of Spifzer and Harm^° in fhe nonrelafivisfic limif. We 
also give analyfical expressions for fhe conducfivily for fhe limiting cases Z ^ oo 
and Z ^ Oin Eqs. (41) and (42). 
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Appendix A. Properties of the homogeneous solutions 

In this appendix we will develop some properties of the solutions of the homo¬ 
geneous radial equations (17). 


1. Definitions 


First we shall obtain fundamental solutions and to the lowest- 

order homogeneous differential equation 


Ll,aX{z) = 0, 


where 


d'^X 


Ll^aXi^) — (1 + ^ \ 


dx ( l{l + 1 ) 

dz \ 


(Al) 


+ 0^-1 X- (A2) 


The variable z eorresponds to rr/c in the main text. In this appendix, I must be an 
integer, a must be real, and z is in the eomplex plane eut along the negative real 
axis. 

By ehanging the independent variable to 7 = \/l + z^ and defining a new 
dependent variable ^( 7 ) = y/zx{z), we obtain the equation 




/'2 10 d^ 

(7 -1)— + “ 

dx 


(a + i)(a - 4) + 




72 _ 1 


e = o. 


This is the generalized Legendre equation, whose solutions are the assoeiated Leg¬ 
endre funetions and choose to define and 

by 

7 qi]a(^) = (A3) 

vimJM = (A4) 

An alfernafive represenfafion is obfained by Whipple’s fransformafion (see Eq. 
(8.739) of Gradshfeyn and Ryzhik^) 


1 


Qiix/z), 


ji[i]aiz) zY{a + l + l) 
1 t_iyp-ta7r 
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The Wronskian of the pair and follows from Eq. (8.741) of Gradshteyn 
and Ryzhik^; for integer I and arbitrary a it is given by 

d d I 

dl[^]a yi[l]a 

Therefore, the solutions and are independent. When I is a non-negative 
integer, then is regular at the origin while yqi]* is singular. 

In order to express the solutions to the higher-order equations, we introduee 
the funetions and yi\k]*{z) where the subseript stands for a string of k 

indiees. The funetions are defined reeursively by 


/ , _ JZ[fc-|-l]*a('^) 7qA:+l]*a'('^) _ 

3l\k-\-2Waa'\^) o To 5 

the case a = a' is handled by taking the limit 

^ _ ^J/[/c+l]*a ('^) 

Jl\k+2]*aa\^) K7 ■ 


(A6a) 


(A6b) 


The funetions yi[k]*{z) are defined in fhe same way, replaeing j by y fhroughouf. 
Making use of fhe properfy {Li^a — Li,a')x = if is readily esfablished 

fhaf and yi[k]^ salisfy 


Ll,ayi[k+l]*a{z) = yi[k]*{z)- (A7b) 


The funetions ji[k]* and yi[k]* are invarianf under permufafion of fhe indiees in 
* and invarianf under ehange of sign of any index in fhese properfies refieef fhe 
eommufivify of L; ^ and L; ,j/ and fhe symmefry L; „ = L; Also, we ean write 

yi[k]* = (A 8 ) 

whieh relleels Li^a = L-i-i,a, and finally we have ,,(- 2 ;) = (-l)'y)[fc]*( 2 ;) and 
yi[k]*(.~z) = {—^yyi[k]*{z)- In mosf of whaf follows, we lisf only fhe properties 
of Equation (A8) may be used fo give fhe eorresponding properties of yi[k]*- 


2. Taylor series 

In fhe limif 2 ; ^ 0, we ean expand (t) in a Taylor series using formula 

(3.2.20) of Bafeman^^ 

^ — r(i -^) —+ yA- 7^), 
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where F is the hypergeometrie funetion 


^ ab z a(a + l)b(b + 1) z'^ 

F{a,b-c-,z) = l + --+ ^ ^ \ - — + ■■■■ 

c 1! c(c + 1) 2! 

Substituting this into Eq. (A3), we obtain a Taylor-series expansion of ji[i]a{z 

^( + 1 —a( + l + a, 3 ^ 


Jl[l]c 


[Z = 


( 2 / + 1 )!! 


F 




(A9) 


where n!! denotes a double faetorial 


n!! = 


2n/2p( 

^n + 1), forn non-negative and even, 
+ 1), forn odd. 


With this definition, we find (2n — !)!!(—2n — 1)!! = (—1)”. 

Anofher useful Taylor series is obfained by a fransformafion of fhe hypergeo- 
mefrie funelion (Eq. (9.131) of Gradshfeyn and Ryzhik^) 


F{a,b-,c;z) = {l — zy “ ^F{c — a,c — b;c-, z). 


This gives 

Jz[i]a(^) _ z^ /f + 2 — af + 2 + a 3 2^ 

7 ~ J¥TT)V.^ \ 2 ’ 2 ’ )' 


(AlO) 


Series (A9) and (AlO) eifher ferminafe (i.e., eonverge everywhere) or eonverge 
inside fhe eirele \z\ = 1. In parfieular, if a is an infeger and I < |a|, fhen one of fhe 
series will terminate. 

It is not so easy to write down Taylor series for the higher order funetions 
ji[k]*{z) (k > 1). However, the leading term is independent of the indiees making 
up *, and ean be obtained by substituting a series solution into Eqs. (A7). One finds 


jl[k]*{z') 

yi[k]*iz') 


Z 


l+2k-2 


(2A:-2)!!(2f + 2A:-l)!! ^ 

(-l)^(2f-2fc + l)!! 

(2A:-2)!!2'-2fc+3 k 


(Alla) 

(Allb) 
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3. Asymptotic series 

Since depends only on the magnitude of a, we can, without loss of gen¬ 
erality, take a > 0 in this section. In the limit z ^ -|-oo, formula (3.2.21) of 
Bateman^^ may be used to derive the following asymptotic expression for 


il|l|a (22)“ 'r(( + a+l)^( 


I — a + 1 —I — a 
2 ’ 2 ’ 

f / “h n -f 1 — I CL 


—a -|- 1; —K 


;a + l;-^ .(A12) 


This cannot be used directly when a is an integer or / ± a is a negative integer; 
in those cases, a limit must be taken. For a and I both non-negative integers, the 
leading order behavior is found to be 


(a-1)! 


{2z 


\a—l 


Jl[l]c 


(a + l)l 

sinh“^ z — X]i=i 
iTz 


for a / 0, 
, for a = 0, 


(A13) 


and 


yi[i]c 




2(a + 0! 


(A14) 


a! 


{2z 


a—1 


for I > a. 


For 0 < I < a and a an integer, yi[^a{z) have the same behavior as 

z ^ -|-(X) (namely z°‘~^). It is, therefore, convenient to define some combination 
of 3i[i]a and 2/z[i]a in which the leading-order behavior cancels. Such a combination 
is 

C \ ■ C \ C 1 ~ ^ ~ f )• ( \ 

- (-1) (a + l)\ 




(A15) 


For I < —a, we have qi\^^a = jz[i]a; for I > a, qi[i]a diverges. The asymptotic 
behavior of qi^a for I < a is 


Qi[i]c 


{-iy+^2{a - I - 1)1 


a: 


(2z 


a—1 


(A16) 
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An asymptotic series for ji[i\a{z) valid for I +oo follows from formula 
(3.2.14) of Bateman^^ 

1 / 2 
The leading term is 

Jdi]a(^) ^ Y{f^/2) ^ for I +00- (AIV) 

4. Recurrence relations 

Recurrence relations that give in terms of y)±i[i]a and y)[i]a±i may be 
derived from the corresponding relations for Legendre functions^’ 

z r 

JZ[l]a(^) = ^21 _|_ 2'«-l[l]a(^) 

+ (( + 1 — a){l + 1 + a)ji^i[i]aiz) (A18a) 

= {I — a + l)jz[i]a-i(^) — {I + a + l)y)[i]a+i(^^) (A18b) 

= ^jz-i[i]a(^) -(/ + ! + a)jiii]a+i{z) ■ (A18c) 

These relations may be combined to express ji[i]a in terms of any pair of its neigh¬ 

bors. The derivative of ji[i]a may be found from 

-^jl[l]aiz) = ^jl-l[l]a{z) - ^-^jl[l]a{z). (A19) 

By combining Eqs. (A6) with the relations (A18a) and (A19), we can generalize 
these recurrence relations to multiple indices 

JZ[fe]*(^) JZ—2[fc-|-l]*a('^) (2^ f) ~JZ—l[fe-|-l]*a(^) 

+ — a^)ji[k+i]*a{z)., (A20) 

= ^ji-i[k]*{z) - ^-^^ji[k]*{z). (A21) 


1 1,3 ^2 

— Ci — 1 Qj — A i- 

2 ’ 2 ’ 2 ’ 7 + 1 
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Other recurrences involving the higher order functions may be found by differen¬ 
tiating Eq. (A18c) with respect to a: 

2{a -h !)(/ -h 1 -h a)jz[ 2 ]a-|-l,a-|-l( 2 :) 

= 2a[zji_i[2]a,a{z) “ “ a)-fjl[2]a,a{^)] 

+ ljl[l]aiz) - jl[l]a+liz), (A22a) 

4(a + 1)^(^ + 1 + a)2z[3]a-|-l,a+l,a-|-l(^) 

= 4o^[^2Z-l[3]a,a,a(^;) “ “ a)7Jz[3]a,a,a(^:)] + ^JZ-l[2]a,a(^) 

-{I- 3a)7jq2]a,a(^) - (/ + 3a + 3)jq2]a+i,a-ri(^)- (A22b) 


Equations (A18) and (A 19) constitute a set of second-order recurrence relations 
which ji[i]a solve. There is another independent solution which we will write as 
If we substitute 

5'Z[l]a(^) = Cl^ayi[l]a{z) 

into the recurrence relations, the resulting relations are equivalent to the original 
set provided ci^a satisfies 

Q,a _ 1 

Q_i,a {l + a){l-a)' 


Cl,a 


/“hi — Qj 


Q,a—1 I -\- Qj 

Depending on the value of I, ci^a may be written in one of three forms 

(-ir 


{I - a)\{l + a)\’ 

ci,a=<^ _ (-l)'(a-/-l)! 

Cla — - 


for I > a, 


for —a<l<a. 


{l + a)\ 

I cfa = ^ “ 1)K“® “ ^ “ 1)^ for ^ 


(A23) 


where we have taken a > 0. Thus = cf^yi[i]a{z) solves the recurrence 

relations for — a < / < a. If this solution is extended beyond this range, it degen¬ 
erates to 0 (for I < —a) or infinity (for I > a). The independence of ji[i]a and gi[i]a 
may be verified by computing the Casoratians. These may be found by substituting 
recurrence relations (A18) and (A19) into the Wronskian, Eq. (A5) to give 

9l[l]aiz)jl—i^l^aiz) 2z[l]a(■^)9Z—l[l]a(■^) Cl,a ^2 ’ 


9Z[l]a(^)jZ[l]a—1 (■^) jl[l]a(.z)gi^ija—l{z') /1 a Z 
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The function qi[^a{z) defined in Eq. (A15) may be written as ~ 

g\yi\a^z). Thus qi[i]a{z) solves the recurrence relations for I < a, and this gives 
a solution independent of j/[i]a for — a < t < a. 

5. Generating function 

From the differential recurrence, Eq. (A21), we obtain 

fl+m[k]*i'y) ( 1 ) (^rym 

where 

Performing a Taylor-series expansion of /o[fe ]*(7 — ez) gives 

/o[fc]*(7 -ez) = J2 fi[k]*il)^j^- 
1=0 

This leads to the following generating function 

_ OO I 

yo[k\*{\lz‘^{l + e^) - 2€Z'y^ = XI (A24) 

z=o 


yi[k]*(V7^ - 1 ) 
(7^ — 1)^/^ 


6. Special cases 


Simple closed expressions exist in certain special cases (in the following equa¬ 
tions, cj = sinh“^ z): 


Jalljaiz) 
Ja—l[l]a(^) 
io[l]a(^) 
j-l[l]aiz) 
J—a—l[l]a(^) 
JO[A:]0...o(^) 


1 fz ^'2a 

_ - _ / —dz' 

( 2 a - l)!b“+i 7o 7 ' ’ 

za-i 

(2a-1)!!’ 

sinh(acr) _ (7 + 2 :)“ — (7 + 2 )““ 

az 2az 

cosh(acr) (7 + 2 )“ + (7 + 2 )““ 
2 : 2z 

( 2 a^ 

^2k-l 

{2k-iy.z’ 


(A25a) 

(A25b) 

(A25c) 

(A25d) 

(A25e) 

(A25f) 
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where we have taken a > 0. 

The reeursion relations, Eqs. (A 18), (A20), and (A22), together with Eq. 
(A25f) allow for 0 < A: < 3 for integer I and integer indiees * to be expressed 
in terms of elementary funetions. The multiple-index homogeneous solutions that 
we need ean be expressed simply in terms of the single-index solutions: 


Jl[2]02(^) — ■XJ/+1[1]1) 


(A26a) 


jl[2]ll{z) — -xii+llllO) 


(A26b) 


Jl[2]22[^) — ■XJ/+1[1]1 + ■X‘J/+2[1]0) 


(A26e) 


JZ[3]022(2:) — -^iz+2[l]0- 


(A26d) 


Eor definiteness, we eatalog all the required funetions ji[k]*{z) and yi[k]*{z) for 
I = 0 and 1: 

<7 1 

■?o[i]o = —) yo[i]o = — {Alla) 

iofili = 1) yofili = (A27b) 


Jo[i ]2 — 7) yo[i ]2 


1 -h 2^2 


(A27e) 


70[2]02 


70[2]11 


70[2]22 


z'^ — a 
4z ’ 

7 CJ — z 
2 z ’ 

—z'y -h (1 -|- 2z‘^)a 

si ’ 


J0[3]022 


—3^7 -|- (3 -|- 2z‘^)a 
32i ’ 


2/0 [ 2 ] 02 


2/0 [ 2 ] 11 


2/0 [ 2 ] 22 


2/0 [3] 022 


— 7 CJ + z 


(A27d) 


(A27e) 


(A27f) 


(A27g) 


'ya — z 

2 ’ 
Z^ 

z^ — a 


( 1 -2^2)7 


(A28a) 


(A28b) 


(A28e) 
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Jl[2]02 = 

-37cr -h {3z -)- z-^) 

7 

yi[2]02 — 

(A28d) 

12^2 

Jl[2]ll = 

—3z'y -|- (3 -|- 22 :^)cj 

1 

yi[2]ii — 2’ 

(A28e) 

8^2 

Jl[2]22 = 

— (3 — 62 ;^) 7 cj -|- {3z — 5z^) 

'-) + za 

(A28f) 

72^2 

yi[2]22 — 2 ’ 

it [3] 022 = 

(15 -|- 62 ;^) 7 cj — (ISz -|- llz^) 

za 

(A28g) 

288^2 

yi[3]022 = 


The function yo[fc]*(' 2 ) is the kernel that occurs in Eq. (10). 

7. Numerical methods 

Our final concern is to determine a method by which ji[k]*, yi[k]*’ their 
derivatives may be calculated accurately and quickly. A direct use of the analytic 
forms is ill-advised, particularly when z is small and I is large. For example, con¬ 
sider the numerator of the analytic form for yi[ 2]02 given in Eq. (A28d) in the limit 
z ^ 0. This consists of three terms, the largest of which is proportional to z; 
however, the sum is proportional to z^. 

Given y)[i]a, Eqs. (A26) may be used to calculate the required multiple-index 
solutions. Furthermore Eq. (A21) may be used to give the derivatives of ji[k]*- 
Equation (A4) gives in terms of ji[k]*- Thus the problem is reduced to 
calculating of all integer —L — 1 < ( < L, integer a > 0 , and real 

z > 0. (For a < 0, we can use ji[\]-a{z) = ji[i]a{z)- For z < 0, we can use 
= (-l)^2)[l]a(^)-) 

We will also be able to avoid problems with numerical underflow and overflow 
by computing ^ where 


Z'' 

ji[i]aiz') (2( -|- jn' 

In the nonrelativistic limit z ^ 0, we have, from Eq. (A9), ^ ^ 1. 

Our main tool for calculating y)[i]a will be the recurrence relation (A18a) which 
when written in terms of j; „ becomes 


Jl-2,a = — 


{I — a){l + a) 

{21 - 1)(2/ -h 1) 




(A29) 
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In order to apply it we need to examine the stability of the reeurrenee for large 1. 
In this limit, the reeurrenee relation is approximately 

Jl—2,a ~ 4 ^ Jl,a: 


whose solution is 

Jl,a ^ 2 

Jl-l^a 7 ± 1 ’ 

Comparing this with the leading term in the asymptotie series (A 17) shows that the 
solution we want corresponds to the upper sign. This solution is dominant when 
the recurrence relation is applied in the backwards direction. If we start with large 
I with some arbitrary mixture of the dominant and subdominant solutions, then on 
each application of the backwards recurrence relation, the subdominant solution 
decreases by (7 — l )/(7 + l) = 2;^/(7 + 1 )^ relative to the desired solution. 

If we desire to compute j/^ to accuracy (5 for / < L, we choose an L' such that 


L' > L + 


21 og( 2:/(7 + 1 ))' 


We set and „ so that their ratio is given by = 2/(7 + 1 )- 

We then use Eq. (A29) as a backwards recurrence to give „ for 0 < / < L. At 
this point j; ^ differs from the desired solution only by an overall multiplicative 
factor. This may be determined from Eq. (A25b) which gives ja-i,a = 1- Thus 
Ji,a = j'l a/ja -1 a- Bccausc of the degeneracy in Eq. (A29) the values of ji^a for 
0 < ^ < a are independent of the choice of L' and starting values ^ and J //_4 
The recurrence is effectively restarted at ^ = a — 1. 

Various optimizations to this scheme are possible. Eor example, it is only nec¬ 
essary to start the recursion at I = L' for one value of a, e.g., a = 0. Eor other 
values of a we can start the recursion at / = L by rewriting Eq. (A 18c) as 


3l,a 


(2/ + l)j/_i^a_l — (/ + 1 

I CL 


a)l3i,a-i 


(A30) 


and using this recurrence to give jL,a and jL-i,a in terms of jL,a-i, 3L-i,a-i, and 
3L-2,a-l- 

Eor large values of z, L' becomes large because the behavior of the dominant 
and subdominant solutions is nearly the same. It is, therefore, possible to use 
forward recursion using Eq. (A29) to obtain j/„ for I > a. Starting values are given 
by Ja-i,a = 1 and ja,a which may be calculated using Jo,o = ojz and recurrence 
relation (A30). As before, backwards recursion should be used for 0 < / < a. 
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For I < —a, we can compute j; „ by backwards recursion using Eq. (A29) 
together with the starting value = 1- For —a < / < 0, we could con¬ 

tinue the backwards recursion using as starting values Jo,a and Ji,^ as found above. 
However, it is sometimes useful to be able to compute qi[i]a', but this cannot be 
accurately computed using Eq. (A 15) when z is large because of the large can¬ 
cellation that occurs in this limit. Instead, we compute (?/[!]„ directly. To avoid 
problems with underflow and overflow, we work with qi^a which is defined by 

Qi[i]aiz) = 

We have seen that qqija satisfies the same recurrence relations as Jqija- This implies 
that qi^a satisfies the recurrence relation (A29). Erom Eqs. (A25c) and (A25d), 
together with the recurrence relation, we have 


=h + z)<.- 

(A3 la) 

+ + 

II 

1 

(A3 lb) 


We can then utilize backwards recursion using Eq. (A29) to give for all —a < 
I < 0. Einally, we can compute J;,^ for —a < I < 0 using Eq. (A15) which gives 




(A32) 
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Table 


TABLE I. Conductivities for various values of the normalized temperature 0 = 
Te/(511 keV) and the effective ion charge state Z. The conductivities are normal¬ 
ized to log A®/®)(Te'^^/Z). 


\ z 
0 \ 


1 

2 

5 

10 

oo 

0 

3.75994 

7.42898 

8.75460 

10.39122 

11.33006 

12.76615 


3.75490 

7.27359 

8.53281 

10.07781 

10.95869 

12.29716 


3.74920 

7.12772 

8.32655 

9.78962 

10.61952 

11.87371 

0.05 

3.72852 

6.73805 

7.78445 

9.04621 

9.75405 

10.81201 

0.1 

3.68420 

6.20946 

7.06892 

8.09361 

8.66306 

9.50746 

0.2 

3.57129 

5.43667 

6.06243 

6.80431 

7.21564 

7.82693 

0.5 

3.18206 

4.13733 

4.47244 

4.88050 

5.11377 

5.47602 

1 

2.65006 

3.13472 

3.32611 

3.57303 

3.72206 

3.96944 

2 

2.03127 

2.27862 

2.39205 

2.54842 

2.64827 

2.82473 

5 

1.33009 

1.45375 

1.51805 

1.61157 

1.67382 

1.78870 

10 

0.94648 

1.02875 

1.07308 

1.13856 

1.18263 

1.26490 

20 

0.67042 

0.72743 

0.75853 

0.80472 

0.83593 

0.89443 

50 

0.42422 

0.46003 

0.47965 

0.50885 

0.52861 

0.56569 

100 

0.29999 

0.32528 

0.33915 

0.35979 

0.37377 

0.40000 
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